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Abstract 

Einstein's theory of gravitation that governs the geometry of space-time, coupled 
with spectacular advance in cosmological observations, promises to deliver a 'stan- 
Oh' dard model' of cosmology in the near future. However, local geometry of space con- 

strains, but does not dictate the topology of the cosmos, hence. Cosmic topology has 
remained an enigmatic aspect of the 'standard model' of cosmology. Recent advance 
' in the quantity and quality of observations has brought this issue within the realm 

of observational query. The breakdown of statistical homogeneity and isotropy of 
cosmic perturbations is a generic consequence of non trivial cosmic topology arising 
^ ' from to the imposed 'crystallographic' periodicity on the eigenstates of the Laplacian. 

' The sky maps of Cosmic Microwave Background (CMB) anisotropy and polarization 

. most promising observations that would carry signatures of a violation of statisti- 

■ cal isotropy and homogeneity. Hence, a measurable spectroscopy of cosmic topology 

0^ ' is made possible using the Bipolar power spectrum (BiPS) of the temperature and 

, polarization that quantifies violation of statistical isotropy. 

o 

qH, I. INTRODUCTION 

U 

bJO' I feel honored to be invited to contribute to this volume honoring the memory of Professor Amal 

Kumar Raychaudhuri (fondly known as AKR) - a great scientist and teacher. The Raychaudhuri 
equations describing the evolution of anisotropic universe models are the footprints of homegrown Indian 
' science in the field of cosmology. Though the background universe is observationally consistent with 

homogeneous and isotropic Friedmann models, the Raychaudhuri equations appears in the evolution of 
inhomogeneities that led to the formation of large scale structures in the universe. It is fair to say 
much of the recent progress in cosmology has come from the interplay between refinement of the theories 
of structure formation and the improvement of the observations. Hence, the Raychaudhuri equations 
have remained as relevant and ingrained in contemporary cosmology as when first put forward by AKR. 
This article describing our ongoing research determine the topology of the universe from the exquisite 
measurements of anisotropy in the Cosmic Microwave background is a humble tribute to the doyen of 
Indian science. 

The realization that a universe with the same local geometry has many different choices of global 
topology has been a theoretical curiosity as old as modern cosmology. De Sitter was quick to point out 
that the first modern model of the cosmos, Einstein's closed {S^, spherical geometry and static) universe 
model, could equally well correspond to the multiply connected 'Elliptical' universe where antipodal 
point of are topologically identified {S^/Z2)- Fig. H depicts the prevalent modern view within the 
concept of inflation, that this relatively smooth 'Hubble volume' that we observe is perhaps a tiny 
patch of an extremely inhomogeneous and complex spatial manifold. The complexity could involve non- 
trivial topology (multiple connectivity) on these ultra-large scales. Given the observational support for a 
homogeneous Hubble volume around us, the diverse possibility of global structure reduces to the tractable 
study limited to spaces of uniform curvature (locally homogeneous and isotropic FLRW models) but with 
non-trivial cosmic topology. For example, in Fig. ^ observers in the 'handle' regions would perceive an 
open (hyperbolic geometry ) universe, and those in the 'bulb' region observe a closed (spherical geometry) 
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universe. Although, in a generic manifold, sectors with exact Euclidean geometry are rare, an epoch of 
inflation in the early universe inflate any region to a nearly flat, Euclidean geometry. 

The motivation for non-trivial topology and the quest to determine size and the shape of our universe 
has a rich and diverse history in modern cosmology Q, |^ |M ^ |^ . Spatially compact universe models 
have both deep theoretical and philosophical appeal - eg., the quantum creation of finite sized universe 
from vacuum is a theoretical motivation and a philosophical abhorrence of any 'infinity' in nature 
would argue against an infinite sized universe 0. A compact universe, with the exception of the three 
sphere, implies a multiply connected universe (i.e., nontrivial cosmic topology). 

The photons of the cosmic microwave background(CMB) propagate freely over distances comparable 
to the cosmic horizon. The CMB anisotropy and its polarization are the most promising observational 
probes of the global spatial structure of the universe on length scales near to and even somewhat beyond 
the 'horizon' scale cHq^). A generic consequence 136] of cosmic topology is the breaking of statistical 
isotropy in characteristic patterns determined by the photon geodesic structure of the manifold . The 
increasingly exquisite measurements of CMB anisotropy have brought cosmic topology from the realm of 
theoretical possibilit y to within the grasp of observations and has received considerable attention over the 
past few years P, 1^ U^, O, H 13 Q ■ CMB polarization not only augments the CMB temperature 
anisotropy observations but is expected to allow a more incisive study of cosmic topology since it arises 
only at the surface of last scattering (SLS). The recent full sky measurements of CMB polarization maps 
by WMAP has opened the door to this new arena (iTj . 



time 




FIG. 1: A cartoon depicting a prevalent view within the inflationary paradigm. The observable universe corresponds to a 
small patch of a very comphcated manifold that has been blown to cosmological scales during an inflationary epoch. Ultra- 
large scale structure could be observable if the the size of this patch is not much smaller that the scales of inhomogeneity 
and non-trivial topology. 

Any multiply-connected space is equivalent to a simply connected space (universal cover) that is tiled 
by Dirichlet domains (DD) under a free acting subgroup, F of the isotropy group, G" of universal cover. 
Hence, even the most complicated cosmic topology reduces to a study of standard FLRW models with 
periodic boundary conditions on the appropriate DD. A simple example being the 3-torus T"^, that is 
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equivalent to the study of a universe in a box with periodic boundary conditions - a routine approximation 
in numerical simulations of large scale structure in the universe. More formally, the is obtained by 
tiling 3-D Euclidean space under the isometry subgroup of discrete translations in three directions. In 
cosmology, the Dirichlet domain constructed around the observer represents the universe as 'seen' by the 
observer. The SI breakdown is apparent in the principal axes present in the shape of the DD constructed 
with the observer located at the base-point (18^ . Equivalently, the fields defining perturbations are built 
of the subset of eigenstates of Laplacian that invariant under F. The correlations of CMB fluctuations in 
the sky would have patterns that are no longer invariant under rotations. 

In a cosmological model with trivial topology, the CMB anisotropy signal is expected to be statisti- 
cally isotropic, i.e., statistical expectation values of the temperature fluctuations AT{q) are preserved 
under rotations of the sky. In particular, the angular correlation function C{q, q') = {AT{q)AT{q')) 
is rotationally invariant for Gaussian flelds, i.e., C{q, q') = C{q ■ q'). In spherical harmonic space, 
where ATiJj) ~ '^im'^i-'mYim{Q)i the condition of statistical isotropy (SI) translates to a diagonal 
{aimaii„i') — C'ldii'Smm' where C/, the widely used angular power spectrum of CMB anisotropy. The 
correlation patterns in the CMB anisotropy that lead to violation of SI implies that imply (a;ma*„) has 



off-diagonal elements. Figure |21 taken from 
for two compact universe models. 



|13| shows the off-diagonal elements in the CMB correlation 
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FIG. 2: The figure taken from 13] illustrates the non-diagonal nature of the expectation values of pair products when 
the CMB anisotropy violates SI in two model compact universe. The radical violation in the model on the left corresponds 
to a small compact universe where CMB photons have traversed across multiple times. The model on the left with mild 
violation of SI corresponds to a universe of size comparable to the observable horizon. For more details, see Il3l 



The observed CMB sky AT(n) is a single realization of the underlying correlation, hence the detection 
of SI violation or correlation patterns pose a great observational challenge [s^l . For statistically isotropic 
CMB sky, the correlation function 

C(ni,n2) EE C(ni • fia) = ^ J dU C(7^nl, Tefia), (1) 

where TZfi denotes the direction obtained under the action of a rotation TZ on n, and dJZ is a volume 
element of the three-dimensional rotation group. The invariance of the underlying statistics under rotation 
allows the estimation of C{hi ■ 112) using the average of the temperature product AT{fi)A.T{n') between 
all pairs of pixels with the angular separation 9. In the absence of statistical isotropy, C(n, n') is estimated 
by a single product A.T{fi)A.T{fi') and, hence, is poorly determined from a single realization. 

Although it is not possible to estimate each element of the full correlation function C{n,n'), some 
measures of statistical anisotropy of the CMB map can be estimated through suitably weighted angular 
averages of AT(h)AT{n'). The angular averaging procedure should be such that the measure involves 
averaging over sufflcient number of independent 'measurements', but should ensure that the averaging 
does not erase all the signature of statistical anisotropy. Recently, we proposed the Bipolar Power 
spectrum (BiPS) Ki {£ = 1,2,3, . . .) of the CMB map as a statistical tool of detecting and measuring 
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departure from SI 19] . The BiPS is formally defined as 

= {2l + lf j dQn, J diln, j dnxH'R)C{nhi, 7^n2)]'. (2) 

In the above expression, C{TZni, TZn2) is the two point correlation at TZfii and TZh2 which are the 
coordinates of the two pixels ni and 712 after rotating the coordinate system by element TZ of the rotation 
group. 

X^{T^) is the trace of the finite rotation matrix in the £M-representation 

I 

x\n) = DiiMm, (3) 

which is called the characteristic function, or the character of the irreducible representation of rank i. It 
is invariant under rotations of the coordinate systems, in eq.JSJ dTZ is the volume element of the three- 
dimensional rotation group. For a statistically isotropic model C(ni, 712) is invariant under rotation, and 
therefore C{TZni, TZn2) — C(ni, 112) and the orthonormality of x^(a;), we will recover the condition for 
SI, 

= K°Sm. (4) 

The Bipolar power spectrum gets it name from its interpretation in the harmonic space. The two 
point correlation of CMB anisotropies, C(ni, 77.2), is a two point function on S'^ x 5^, and hence can be 
expanded as 

C(ni,n2)= {YiAni)®Yi,{n2)}iM, (5) 

where Af^f^ are coefficients of the expansion (here after BipoSH coefficients) and {Yi-^{ni) ® YiaC^-a)}^*/ 
are the Bipolar spherical harmonics which transform as a spherical harmonic with £, M with respect to 
rotations '2(| given by 

{Yi,in,)®Yi,in2)}eM ^ E ^"™i;2m.^i™i("2)5^.m.(^2), (6) 

mi m2 

in which Cf*^^;^^^ are Clebsch-Gordan coefficients. We can inverse-transform C(ni, 712) to get the Af^*^ 
by multiplying both sides of eq.© by {Yi'^{ni) (g) (7T-2)}|'m' ^.nd integrating over all angles, then the 
orthonormality of bipolar harmonics implies that 

4t = / rff^ni / dn^, Cin,, h2) {Yi,{ni) ® >i,(n2)},*M- (7) 

The above expression and the fact that C{fii, 712) is symmetric under the exchange of fii and 712 lead to 
the following symmetries of A\-^f^ 



= (-l)('^+'-^)Af*^„ (8) 

IM X 



^ A\tU,,2U, fc^O, 1,2, ••• 



The Bipolar Spherical Harmonic (BipoSH) coefficients, ^fjj^' ^'"^ linear combinations of off-diagonal 
elements of the harmonic space covariance matrix, 

mim2 

This means that A^^^^ completely represent the information of the covariance matrix in harmonic space 
{o-iimiOLi^m^) ■ When SI holds, the harmonic space covariance matrix is diagonal and hence 

A{{(' = {-l)'Ci{2l + lfl^6u'5tohm, (10) 
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BipoSH expansion is the most general representation of the two point correlation functions of CMB 
anisotropy. The well known angular power spectrum, Ci is a subspace of BipoSH coefficients correspond- 
ing to the that represent the statistically isotropic part of a general correlation function. When SI 
holds, or equivalently C; have all the information of the field. But when SI breaks down, A^^ are not 
adequate for describing the field, and one needs to take the other terms into account. The Bipolar power 
spectrum (BiPS) is defined as a rotationally invariant contraction of the BipoSH coefficients 



This definition is identical to the real space expression in eq. I^Jl. More importantly, BiPS is mea- 
surable from a single CMB map since averages over many independent modes and reduces the cosmic 
variance |3i | 1,19]. The BiPS of the CMB anisotropy maps measured by WMAP has been computed |2l|. 
Preliminary BiPS results on the CMB polarization maps from the three year of WMAP data have also 
emerged in past few months [23, 113 ■ 

The BiPS is sensitive to structures and patterns in the underlying total two-point correlation function 
[l9ll2^ . The BiPS is particularly sensitive to real space correlation patterns (preferred directions, etc.) 
on characteristic angular scales. In harmonic space, the BiPS at multipole t sums power in off-diagonal 
elements of the covariance matrix, {aimai'm'), in the same way that the 'angular momentum' addition of 
states Im, I'm' have non-zero overlap with a state with angular momentum |Z — Z'| < i < l + l' . Signatures, 
like aim and ai+nm being correlated over a significant range / are ideal targets for BiPS. These are typical 
of SI violation due to cosmic topology and the predicted BiPS in these models have a strong spectral 
signature in the bipolar multipole £ space The orientation independence of BiPS is an advantage 

since one can obtain constraints on cosmic topology that do not depend on the unknown (but specific) 
orientation of the pattern {e.g., preferred directions of DD relative to the sky). 

Spaces of constant curvature have been completely classified . For Euclidean geometry, there 

are known to be six possible topologies that lead to orientable spaces. The simple flat torus, A4 = , is 
obtained by identifying the universal cover Ai" = £^ under a discrete group of translations along three 
non-degenerate axes, Si, S2, S3: ^ + nLi, where Li is the identification length of the torus along Si 
and n is a vector with integer components. In the most general form, the fundamental domain (FD) is 
a parallelepiped defined by three sides Li and the three angles ai between the axes ('squeezed torus'). If 
Si are orthogonal then one gets cuboid FD, which for equal Li reduces to the cubic torus. The cuboid 
and squeezed spaces which can be obtained by a linear coordinate transformation C on cubic torus can 
have distinctly different global symmetry ^|. 

Study of the BiPS signature of cosmic topology has already been undertaken and ongoing p^l2^ l2^. 
The correlation function C((j, q') for CMB anisotropy in a multiply-connected universe such as the torus 
space can be computed using the regularized method of linages [l^. Fig.^lplots the predicted Ki spectrum 
for a number of cubic, cuboid and squeezed torus spaces [ig- Similar results for Poincare Dodecahedron 
show a characteristic BiPS with a dominant peak at kq over a range of value of curvature radius (including 
integrated Sachs- Wolfe effect). This relates to the angular separation of the directions to faces of the 
extremely symmetric DD of Poincare dodecahedral space. Fig. O shows that K£ is zero for odd £. This is 
intimately related to the symmetries of the Dirichlet domain, which in turn is dictated by the properties 
of the subgroup of isometrics F. The BiPS is zero at odd bipolar multipoles if DD has 2n-fold symmetry 
about an axis and reflection symmetry in an orthogonal plane. It can be proved that all Euclidean and 
all spherical spaces generated by single action F satisfy this condition. Remarkably enough, compact 
hyperbolic spaces do not satisfy these condition, and are generically expected to have non-zero BiPS at 
odd value bipolar multipoles £ |[29|. This provides a measurable classification of cosmic topology based 
on CMB anisotropy and polarization, i.e., a spectroscopy of cosmic topology. 

A simple working example is the BiPS signature of a non-trivial topology can be given for a universe, 
where the correlation function is given by 




(11) 



1,1', M 




(12) 



n 



in which, n is 3-tuple of integers (in order to avoid confusion, we use q to represent the direction instead 
of n), the small parameter < 1 is the physical distance to the SLS along q in units of L/2 (more 
generally, L/2 where L = (LiL2-^3)^^'^) and L is the size of the Dirichlet domain (DD). When e is a small 
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constant, the leading order terms in the correlation function cq. (|12|l can be readily obtained in power 
series expansion in powers of e. For the lowest wave numbers |np = 1 in a cuboid torus 



C{q,q') « 2^P$(2Vi,)cos(7re/3,Ag,) 



(13) 



Cn 



where Aqi are the components of A.q = q — q' along the three axes of the torus and Pi = L/Li. From 
this, the non-zero can be analytically computed to be 



^2(l_4e2 

122887r2 „ 
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(14) 



K4 has the information of the relative size of the Dirichlet domain and one can use it to constrain the 
topology of the universe. 
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FIG. 3: The figure taken from llSl shows the BiPS ki spectra for flat tori models . The top row panels are for cubic tori 
spaces. The left panel shows spaces of volume, Vm ; larger than the volume V* contained in the sphere of last scattering (SLS) 
with = 3.7,1.9,1.1, respectively. The right panel shows small spaces with VmIV* = 0.24,0.07,0.03, respectively. 

Note that K2 = for cubic tori. The middle panels consider cuboid tori with 1 ; 5 and 1 : 2 ratio of identification lengths. 
The bottom panels show Ki for equal-sided squeezed tori with a = 45° , 60° and 75° . In the middle and bottom rows, the 
right panels show the case when radius of SLS, i?, = i?< the in-radius of the space. Here, the SLS just touches its nearest 
images which is at the threshold where CMB anisotropy is multiply imaged for larger ij*. The cases in the left panels of 
lower two rows have V^/V* = 1 and are at the divide between large and small spaces. 

The results of WMAP are a milestone in CMB anisotropy measurements since it combines high angular 
resolution, high sensitivity, with 'full' sky coverage allowed by a space mission. The Wilkinson Microwave 
Anisotropy Probe ( WMAP) observations are consistent with the predictions of the concordance ACDM 
model with scale-invariant and adiabatic fluctuations which have been generated during the inflationary 
epoch [sol f3ll [3^. After the first year of WMAP data, the SI of the CMB anisotropy {i.e. rotational 
invariance of n-point correlations) has attracted considerable attention. Tantalizing evidence of SI break- 
down (albeit, in very different guises) that mounted in the WMAP first year sky maps, using a variety of 
different statistics are expected to persist in the three persist in three year data (see 0| for discussion 
and references.). 

The CMB anisotropy map based on the WMAP data are ideal for testing for statistical isotropy. 
Preferred directions and statistically anisotropic CMB anisotropy have been discussed in literature earlier 
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|33l |3J| . A number of direct searches for signature of cosmic topology have been proposed and carried 
out on early CMB data from COBE-DMR. Full Bayesian likelihood comparison to the data of specific 
cosmic topology models is another approach that has applied to COBE-DMR data 0, 0, ^3- The 
generic features of Kg spectrum are related to the symmetries of correlation pattern. For cosmic topology, 
K{ are sensitive to SI violation even when CMB is not multiply imaged. The orientation independence 
of BiPS is an advantage for constraining patterns (preferred directions) with unspecified orientation in 
the CMB sky such as that arising due to cosmic topology or, anisotro pic cosmology [ssj . Extension of 
BiPS analysis to CMB polarization maps has been studied recently [23, 123 adds a new dimension to the 
spectroscopy of cosmic topology. 

In summary, there are strong theoretical and philosophical motivations for a non-trivial cosmic topology. 
The breakdown of statistical homogeneity and isotropy of cosmic perturbations is a generic feature of 
non trivial cosmic topology. A promising observational approach is to hunt for SI violation in the CMB 
anisotropy. The underlying correlation patterns in the CMB anisotropy and polarization in a multiply 
connected universe is related to the symmetry of the Dirichlet domain. BiPS has the advantage of being 
independent of the overall orientation of the Dirichlet domain with respect to the sky. The pattern of 
SI violation of a cosmic topology leads to a measurable, characteristic Bipolar power spectrum related to 
the principle directions in the Dirichlet domain and symmetries of the two point correlation function. 
The Bipolar power spectroscopy of cosmic topology presents itself as promising pursuit for current and 
upcoming measurements of CMB anisotropy and polarization. 

The author acknowledges years of very fruitful collaboration with Dick Bond and Dmitry Pogosyan. 
Amir Hajian devoted most of his doctoral thesis to this effort and the results mentioned here were jointly 
obtained with him The article also includes recent and ongoing work done with Himan Mukhopadhyay 
and Soumen Basak. 
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